function [p]= satpos(w,deltatime,r,center)
    p   = repmat( center,1,length(deltatime))+r*[cos(w*deltatime)',sin(w*deltatime)',zeros(length(deltatime),1)]';
